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ABSTRACT:  To  better  understand  population  connectivity  of  the  blue  crab  CaUinectes  sapidus 
and  how  it  may  have  been  affected  by  the  Deepwater  Horizon  oil  spill,  we  simulated  larval  disper¬ 
sal  with  a  biophysical  model  of  the  coastal  waters  from  western  Louisiana  to  the  Florida  panhan¬ 
dle.  We  investigated  connectivity  patterns,  intra-annual  variability,  and  potential  impacts  of  the 
Deepwater  Horizon  oil  spill  during  the  spring  and  summer  of  2010.  Overall,  we  found  that  the  Mis¬ 
sissippi  River  delta  (MRD)  is  a  barrier  to  dispersal,  and  that  local  retention  was  high;  of  the  7.7  % 
of  larvae  that  successfully  settled,  37.5  %  returned  to  their  natal  estuary  and  28.5  %  to  an  adjacent 
one.  We  used  network  metrics  to  assess  the  overall  diversity  of  population  connectivity  and  the 
importance  of  individual  estuaries  to  maintaining  connectivity.  The  proportion  of  larvae  that  suc¬ 
cessfully  settle  does  not  significantly  change  during  the  spawning  season,  but  connectivity  among 
estuaries  significantly  declines.  Estuaries  near  the  MRD  were  most  important  for  maintaining  con¬ 
nectivity,  likely  because  they  were  the  primary  source  of  the  few  larvae  that  crossed  the  MRD. 
These  patterns  influence  the  distribution  of  settlement  locations  for  larvae  that  were  potentially 
exposed  to  oil.  A  total  of  38. 1  %  of  the  simulated  larvae  were  potentially  exposed  to  oil,  and  these 
larvae  were  concentrated  on  the  eastern  side  of  the  MRD.  For  some  spawning  events,  up  to  96.3  % 
of  the  larvae  that  successfully  settled  east  of  the  MRD  were  potentially  exposed  to  oil,  which  may 
have  substantial  implications  for  population  dynamics.  These  results  provide  quantitative  predic¬ 
tions  regarding  blue  crab  connectivity  in  the  northern  Gulf  of  Mexico  that  can  be  corroborated 
with  data.  The  predictions  can  be  applied  for  disaster  management  planning  and  for  management 
of  this  environmentally  and  economically  important  species. 

KEY  WORDS:  Lagrangian  particle  tracking  •  Graph-theoretic  metric  •  Shannon  index  • 
Vertex  degree  •  Betweenness  centrality  •  Potential  connectivity 
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INTRODUCTION 

Connectivity  has  been  identified  as  an  important 
driver  of  population  dynamics  in  theoretical  (Hast¬ 
ings  &  Harrison  1994,  Botsford  et  al.  2001,  Hastings  & 
Botsford  2006,  Bode  et  al.  2008,  Kininmonth  et  al. 
2010b)  and  empirical  studies  (Meynecke  et  al.  2008, 

'Corresponding  author:  btjones@mit.edu 


Rooker  et  al.  2008,  Hinrichsen  et  al.  2011)  of  marine 
populations.  Population  connectivity  may  be  defined 
as  the  successful  exchange  of  individuals  among 
spatially  separated  sub-populations  that  constitute  a 
single  metapopulation  (Cowen  et  al.  2007).  It  in¬ 
cludes  the  rate,  extent,  and  spatial  structure  of  such 
exchange  (Cowen  et  al.  2006),  as  well  as  the  tempo- 
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ral  variance  of  these  components.  High  levels  of  con¬ 
nectivity  may  increase  resilience  to  local  overhar¬ 
vesting  and  local  disturbances  (Strathmann  et  al. 
2002),  adaptability  to  climate  change  through  migra¬ 
tion  (Peterson  et  al.  2002,  Harley  et  al.  2006),  and 
adaptability  to  regional  level  shifts  via  gene  flow  and 
genetic  diversity  (Bohonak  1999,  Garant  et  al.  2007). 
High  connectivity  thus  decreases  the  impact  of  local, 
regional,  and  global  disturbances  on  the  population 
(Cowen  &  Sponaugle  2009). 

For  many  marine  organisms,  dispersal  occurs  pri¬ 
marily  during  a  planktonic  larval  phase  (Kinlan  & 
Gaines  2003),  Studying  movements  of  these  larvae  in 
the  field  is  challenging  because  of  their  small  size 
and  high  natural  mortality  rates.  Direct  estimates  of 
larval  connectivity  can  sometimes  be  obtained  using 
mark-recapture  of  artificially  tagged  larvae  (Thor- 
rold  et  al.  2002)  or  by  analyzing  geochemical  signa¬ 
tures  in  calcified  structures  (Thorrold  et  al.  2007). 
However,  applying  these  methods  is  difficult  or 
impossible  for  some  teixa,  particularly  crustacean  lar¬ 
vae,  which  contain  no  calcified  parts  and  molt  fre¬ 
quently.  Indirect  estimates  of  connectivity  can  be 
obtained  using  genetic  markers,  although  the  costs 
required  for  developing  appropriate  genetic  markers 
are  prohibitively  high  for  most  researchers  at  this 
time  (Hedgecock  et  al.  2007),  Lagrangian  particle 
tracking  is  an  increasingly  common  method  for  esti¬ 
mating  potential  connectivity  that  can  be  applied  to 
a  wide  variety  of  taxa.  This  method  simulates  the 
movement  of  individual  particles  ('larvae')  through 
hydrological  flow  fields.  Apphcations  of  Lagrangian 
particle  tracking  to  larval  dispersal  have  been  facili¬ 
tated  by  advances  in  computing  and  the  availabUity 
of  high-resolution  hydrodyncunic  data  sets.  These 
models  have  been  used  by  many  studies  for  a  variety 
of  apphcations.  Examples  include  estimating  the  pro¬ 
portion  of  larvae  locally  retained  (Cowen  et  al.  2000), 
investigating  how  larval  behavior  influences  disper¬ 
sal  distance  (Paris  et  al.  2007,  North  et  al.  2008),  and 
determining  the  optimal  spatial  design  of  marine 
reserves  (Hinrichsen  et  al.  2009,  Costello  et  al.  2010, 
Kininmonth  et  al.  2011). 

In  this  study,  we  coupled  a  Lagrangian  particle¬ 
tracking  model  to  output  from  a  high-resolution 
hydrodynamic  model  to  estimate  potential  connectiv¬ 
ity  of  the  blue  crab  Callinectes  sapidus  in  the  north¬ 
ern  Gulf  of  Mexico  (nGOM).  The  blue  crab  is  an 
abundant,  estuarine-dependent  species  in  the  west¬ 
ern  Atlantic  and  Gulf  of  Mexico  (GOM),  where  it 
supports  profitable  commercial  and  recreational  fish¬ 
eries  (Guillory  et  al.  2001).  Blue  crab  fisheries  were 
once  considered  inexhaustible  because  of  the  crab’s 


high  fecundity  and  rapid  growth  rates,  but  recent 
declines  in  the  historically  strong  Chesapeake  Bay 
stock  has  highlighted  the  need  for  better  manage¬ 
ment  of  this  species  (Miller  et  al.  2005).  In  the  US 
GOM,  the  stock  fluctuates  from  year  to  year  (Van 
Voorhees  2012),  but  there  has  been  no  apparent 
downward  trend  in  blue  crab  abundance,  and  the 
industry  has  recently  been  certified  as  sustainable  by 
the  Marine  Stewardship  Council  (DeAlteris  et  al. 
2012).  Currently,  the  US  GOM  blue  crab  stock  is 
managed  at  the  state  level  as  5  separate  fisheries 
(Texas,  Louisiana,  Mississippi,  Alabama,  and  Flo¬ 
rida),  but  the  spatial  structure  of  the  blue  crab  popu¬ 
lation  and  degree  of  dispersal  between  the  states  is 
unknown.  Decades  of  research  in  the  areas  sur¬ 
rounding  the  Delaware  and  Chespeake  Bays  demon¬ 
strated  that  blue  crab  larvae  are  most  likely  to  recruit 
to  their  natal  bay  (Epifanio  &  Tilburg  2008),  but  con¬ 
nectivity  in  the  nGOM  may  be  different  because 
nearshore  circulation  in  the  nGOM  is  primarily 
driven  by  wind,  outflow  from  the  Mississippi  River, 
and  intrusions  of  the  Loop  Current  and  its  periodic 
eddies  (Ward  1980),  Some  previous  genetic  studies 
have  found  evidence  for  overall  Gulf-wide  panmixia 
(McMillen-Jackson  &  Bert  2004),  but  others  have 
foimd  evidence  for  significant  structure  over  shorter 
spaticd  and  temporal  scales  (Kordos  &  Burton  1993, 
Darden  2004,  Yednock  &  Neigel  2014).  This  suggests 
that  long-distance  dispersal  occurs  frequently  enough 
to  keep  the  populations  panmictic,  but  also  that  dis¬ 
persal  could  be  limited  at  scales  relevant  to  ecologi¬ 
cal  processes. 

We  ran  simulations  over  the  course  of  April  to 
October  2010  over  a  spatial  domain  of  coastal 
waters  from  western  Louisiana  to  the  Florida  pan¬ 
handle.  The  spatial  extent  and  time  frame  allowed 
us  to  estimate  potential  effects  of  the  Deepwater 
Horizon  (DWH)  oil  spiU.  The  DWH  spill,  the  largest 
accidental  marine  oil  spill  in  US  history  (Graham  et 
al.  2011),  started  when  the  DWH  oil  platform,  situ¬ 
ated  80  km  southeast  of  Venice,  Louisiana,  ex¬ 
ploded  and  sank  on  20  April  2010.  An  estimated 
4.0  million  barrels  of  oil  were  subsequently  released 
from  the  broken  wellhead  at  a  depth  of  1500  m  until 
the  well  was  capped  on  15  July  2010.  Much  of  the 
oil  remained  deep,  but  about  50%  rose  to  the  sea 
surface,  forming  shcks  that  were  visible  in  sateUite 
imagery  (Peterson  et  al.  2012).  We  expect  that  blue 
crab  larvae  came  into  contact  with  the  surface  oil 
because  they  are  surface-dwelling  organisms,  found 
at  depths  of  1-3  m  (Epifanio  &  Tilburg  2008),  and 
the  release  of  oil  occurred  during  their  primary 
spawning  season  (Guillory  et  al.  2001). 
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A  variety  of  methods  are  available  for  extracting 
meaningful  information  from  the  output  of  Lagran- 
gian  particle-tracking  models;  we  present  our  analy¬ 
sis  in  the  context  of  graph  theory.  Graph  theory  is  a 
field  of  mathematics  that  analyzes  the  encoding  and 
properties  of  the  connections  among  a  set  of  discrete 
objects,  and  the  field  offers  many  usable  metrics  for 
analysis  of  Lagrangian  model  output.  Graph-level 
metrics,  such  as  diameter,  degree  distribution,  and 
minimum  path  length,  have  been  used  in  a  small 
number  of  marine  studies  to  characterize  the  struc¬ 
ture  of  networks  (Treml  et  al.  2008,  Kininmonth  et  al. 
2010a.  2011,  Holstein  et  al.  2014),  and  are  widely 
used  for  the  assessment  of  connectivity  patterns  in 
other  fields.  Other  metrics  are  used  to  identify  im¬ 
portant  vertices  within  the  network.  These  include 
basic  metrics  such  as  the  number  of  settlers  arriving 
to  (destination  strength)  or  coming  from  (source 
strength)  a  particular  vertex  (Mitarai  et  al.  2009,  Wat¬ 
son  et  al.  2010),  and  measures  of  centrality  that  con¬ 
sider  the  importance  of  a  vertex  to  the  entire  network 
(e.g.  closeness  centrality,  betweenness  centrality, 
eigenvalue  centrality;  Treml  et  al.  2008,  Kininmonth 
et  al.  2010a,  Watson  et  al.  2011).  Although  many  of 
the  graph  metrics  have  equivalent  operations  in  lin¬ 
ear  algebra  or  other  fields  of  mathematics,  present¬ 
ing  the  results  in  the  context  of  graph  theory  retains 
the  ecological  interpretation  of  our  system  as  a  net¬ 
work  of  interconnected  populations.  In  this  study,  we 
generated  a  small  network  with  14  vertices,  where 
the  vertices  represent  estuaries  in  the  nGOM  and  the 
edges  between  vertices  represent  the  probability  of 
larvae  traveling  between  pairs  of  estuaries.  We  gen¬ 
erate  overall  connectivity  matrices  for  the  network 
from  our  simulations.  To  estimate  the  importance 
of  individual  estuaries  in  maintaining  connectivity 
throughout  the  region,  we  calculate  betweenness 
centrality  (the  proportion  of  the  most  probable  path¬ 
ways  that  pass  through  each  vertex)  for  each  estuary. 
To  quantify  the  overall  level  of  connectivity  in  the 
network,  we  use  a  combination  of  average  vertex 
degree,  average  vertex  source  strength,  and  average 
vertex  Shannon  index.  The  vertex  out-degree  and 
vertex  source  strength  are  measures  of  the  route- 
specific  flux  through  a  vertex  (Rayfield  et  al.  2011), 
and  the  Sharmon  index  measures  the  topological 
diversity  of  the  network  (Eagle  et  al.  2010).  Vertex 
degree  is  the  number  of  possible  destinations  for 
particles  that  originate  at  that  vertex.  Source  strength 
is  the  sum  of  these  outgoing  edges  and  represents 
the  probability  that  a  particle  released  at  that  vertex 
will  successfully  settle.  The  Shannon  index  combines 
these  indices  into  a  single  metric  that  quantifies  both 


the  number  and  evenness  of  the  dispersal  probabili¬ 
ties.  It  was  originally  derived  as  an  entropy  measure 
of  the  maximum  amount  of  information  that  can 
be  transmitted  in  a  finite  length  message  (Shannon 
1948),  and  has  been  applied  in  ecology  as  a  measure 
of  species  diversity  (Good  1953,  Margalef  1958)  and 
community  stabihty  (MacArthur  1955).  We  apply  it 
to  measuring  metapopulation  stability.  In  a  graph- 
theoretic  context,  the  Shannon  index  for  a  vertex 
is  maximized  for  a  fuUy  and  uniformly  connected 
vertex.  A  vertex  with  fewer  connections  or  a  non- 
uniform  distribution  of  larvae  traversing  each  out¬ 
bound  connection  will  have  a  lower  Shannon  index. 

Using  these  techniques,  we  describe  and  quantify 
the  patterns  of  population  connectivity  for  C,  sapidus 
in  the  nGOM  and  investigate  how  those  patterns 
may  have  been  influenced  by  larval  exposure  to  oil 
from  the  DWH  oil  spill.  We  examine  the  locations  of 
settlement  of  larvae  that  were  potentially  exposed  to 
oil.  Using  graph  metrics,  we  quantify  the  connectivity 
amd  the  changes  in  connectivity  over  the  course  of  a 
single  season,  and  the  importance  of  different  re¬ 
gions  to  maintaining  connectivity.  Using  sensitivity 
analyses,  we  demonstrate  that  our  main  conclusions 
are  robust  to  assumptions  made  by  our  model. 

METHODS 

We  modeled  larvae  as  planktonic  particles  and 
simulated  their  movement  using  a  Lagrangian  Indi¬ 
vidual  Based  Model  (IBM),  then  assessed  the  poten¬ 
tial  dispersal  patterns  and  destinations  of  larvae 
whose  trajectories  crossed  surface  oil  slicks  from  the 
DWH  oil  spill. 

Hydrodynamic  model 

The  IBM  is  forced  by  oceanographic  data  from 
the  Northern  Gulf  of  Mexico  Nowcast/Forecast 
System  (NGOMNFS).  NGOMNFS  is  a  free  surface, 
data-assimilating,  hydrodynamic  model  based  upon 
NCOM  (Martin  et  al.  2009)  that  solves  the  primitive 
equations  for  incompressible,  hydrostatic  flow  with 
the  finite  difference  method  (http;//www7320.nrlssc, 
navy.mil/IASNFS\_WWW/NGOMNFS\_WWW/ 
NGOMNFS.html;  Ko  et  al.  2008).  The  inputs  to 
NGOMNFS  include  bathymetry  that  is  based  upon 
the  2-min  resolution  NRL  DBDB2  data  set  and  im¬ 
proved  with  1-min  resolution  NOAA  NGDC  data, 
and  NGOMNFS  assimilates  sea  surface  height  from 
satelhte  altimeters  (GFO-1,  Jason-1,  and  ERS-2)  and 
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sea  surface  temperature  from  MODIS  and  AVHRR 
via  the  NRL  MODAS  (Fox  et  al.  2003).  Three-hourly 
wind  stress,  sea-level  air  pressure  and  heat  fluxes 
from  NOGAPS  and  COAMPS  atmospheric  prediction 
systems  are  applied  for  surface  forcing.  Tidal  forcing 
and  freshwater  discharges  from  112  rivers  or  runoff 
points  are  also  applied,  and  open  boundary  condi¬ 
tions  are  provided  by  the  large-scale  Intra- Americas 
Sea  Nowcast/Forecast  System  (Ko  et  al.  2003). 
NGOMNFS  produces  hourly  1/72“  (-1.9  km)  3D  cur¬ 
rent  velocity,  temperature,  and  salinity  nowcasts  for 
the  nGOM  above  27.5“ N  and  from  92“  W  to  86“  W  at 
37  vertical  sigma  levels.  NGMONFS  output  is  avail¬ 
able  from  2001  until  the  present.  Validation  with 
NOAA  tidal  observations  shows  that  NGOMNFS 
performs  very  well  in  the  prediction  of  tides  and 
subtidal  wind-driven  sea-level  variation  (r  -  0.9),  and 
NGOMNFS  output  has  been  applied  in  various 
studies  (e.g.  D'Sa  &  Ko  2008,  Green  et  al.  2008). 

Individual-based  model 

The  IBM  advects  and  diffuses  particles  in  3  dimen¬ 
sions,  and  implements  swimming  behavior  that 
restricts  larvae  to  the  surface  layer.  Particle  trajecto¬ 
ries  are  adverted  using  a  4th-order  Runga-Kutta 
scheme  and  fixed  hourly  time  step.  Within  each  grid 
cell,  velocities  are  interpolated  to  the  particle  posi¬ 
tion  using  a  tri-linear  interpolation.  Although  linear 
interpolations  are  divergent,  the  NGOMNFS  output 
is  of  high  resolution,  and  the  interpolation  errors  are 
minor  relative  to  the  broader  problem  of  how  weU  the 
ocean  model  simulates  the  ocean  dynamics.  There¬ 
fore,  we  used  the  computationally  efficient  linear 
interpolation  scheme.  Sub-grid-scale  diffusion  is  sim¬ 
ulated  using  a  random  flight  model,  as  has  been  done 
in  similar  studies  (Marinone  et  al.  2004,  Johnson  et  al. 
2013).  The  random  flight  model  assumes  that  sub- 
grid-scale  turbulence  is  isotropic  in  the  horizontal 
direction.  It  uses  normally  distributed  random  sam¬ 
ples,  the  time  step  duration,  and  a  diffusion  coeffi¬ 
cient  to  adjust  the  particle  position  in  the  x  and  y 
directions  at  each  time  step.  We  chose  a  diffusion 
coefficient  of  1.21  m^  s"'  based  upon  the  linear  rela¬ 
tionship  between  length  scale  and  diffusion  pre¬ 
sented  by  Okubo  (1971)  and  the  characteristics  of 
NGOMNFS.  In  addition  to  the  passive  advection 
and  diffusion  of  particles,  we  implemented  vertical 
swimming  behavior  that  maintains  them  near  the 
surface  because  field  and  laboratory  studies  show 
that  the  larvae  inhabit  surface  water  throughout  their 
larval  development  (Epifanio  1988,  Epifanio  &  Til¬ 


burg  2008).  In  the  model,  if  the  depth  of  a  larva  was 
at  least  2  m,  then  a  simple  linear  function  of  z,+i  = 
O.lz,  was  apphed,  where  z,  is  the  vertical  position 
of  the  larva  at  time  t.  This  formulation  assumes  that 
larvae  swim  upward  more  forcefully  in  response  to 
increased  downward  transport.  This  function  is  con¬ 
sistent  with  the  zoeal  swimming  speeds  of  0.005  m  s"' 
measured  by  Sulkin  et  al.  (1980). 

Using  this  IBM  together  with  spawning  and  larval 
behavior  parameters  from  the  known  natural  history 
of  blue  crabs,  we  ran  24  simulations  of  blue  crab 
dispersal.  Exact  spawning  locations  are  mostly  un- 
knovra  in  the  nGOM  (but  see  Gelpi  et  al.  2009),  so  we 
released  larvae  throughout  the  plausible  region  of 
potential  sites.  We  identified  locations  within  the 
NGOMNFS  domain  of  3-10  m  depth  and  at  least 
22  ppt  salinity,  which  are  the  known  conditions  nec¬ 
essary  for  blue  crab  spawning  (Guillory  et  al.  2001, 
Rittschof  et  al.  2010).  The  release  points  were  a 
uniform  grid  of  200  potential  spawning  locations 
throughout  this  region  (Fig.  1).  Assessing  population 
connectivity  patterns  requires  us  to  quantify  both 
larval  dispersal  and  the  movements  of  ovigerous 
females  from  estuaries  to  spawning  habitat.  During 
the  analyses,  we  assume  that  each  spawning  site 
may  contain  spawning  females  from  any  estuary 
within  100  km,  and  that  the  proportion  of  females 
from  each  estuary  is  inversely  proportional  to  the 
distance  between  the  estuary  and  the  spavraing  site. 
For  example,  a  spawning  location  that  was  75  km 
from  estuary  i  and  25  km  from  estuary  j  would  be 
assumed  to  have  25  %  of  its  larvae  sourced  from  i  and 
75%  from  j.  Because  the  choice  of  migration  distance 
for  ovigerous  females  was  arbitrary,  we  performed  a 
sensitivity  analysis  to  this  parameter. 

Each  of  the  24  simulations  began  at  the  local  noc¬ 
turnal  ebbing  tide  (Provenzano  et  al.  1983)  on  the 
date  of  each  quarter  of  the  moon  phase,  to  match  the 
observed  highest  spawning  events  (Ziegler  2002) 
from  1  April  2010  until  30  September  2010,  which  is 
the  primary  spawning  season  for  blue  crabs  in  the 
nGOM  (Guillory  et  al.  2001).  For  each  spawning  site, 
tide  data  were  obtained  for  the  closest  location  avail¬ 
able  on  the  NOAA  Tides  and  Currents  database 
(https://tidesandcurrents.noaa.gov). 

At  each  of  the  200  spawning  sites,  we  simulated  the 
release  and  subsequent  hourly  movements  of  10000 
super-individuals.  Each  super-individual  is  assumed 
to  represent  a  large  number  of  larvae,  which  has 
been  used  in  prior  studies  where  the  number  of  lar¬ 
vae  spawned  each  year  greatly  exceeds  computa¬ 
tional  capacities  (Scheffer  et  al.  1995,  Woods  2005, 
Tian  et  al.  2009).  Larvae,  once  spawned,  are  assumed 
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Fig.  1.  The  study  region.  The  Northern  Gulf  of  Mexico  Nowcast/Forecast  System  model  includes  the  coastal  waters  north  of 
27.5°  N  eind  between  92°  W  and  86°  W.  The  14  labeled  estuaries  represent  our  division  of  adult  habitat.  We  assume  that  oviger- 
ous  female  crabs  {Callinectes  sapidus)  migrate  from  these  estuaries  to  the  spawning  sites  that  are  indicated  by  blue  dots.  The 
boxes  surrounding  each  estuary  depict  the  settlement  area.  Additional  details  about  the  spawning  and  settlement  criteria  are 

provided  in  the  'Individual-based  model'  section  of  the  text 


to  remain  within  the  plankton  for  at  least  27  d,  which 
is  the  average  larval  development  period  observed 
in  the  laboratory  (Zmora  et  al.  2005).  Larvae  are 
deemed  to  settle  if  between  27  and  60  d  after  spawn¬ 
ing  they  are  in  water  with  a  depth  of  less  than  0,5  m. 
We  assume  that  all  post-larvae  (megalopae)  in  shal¬ 
low,  nearshore  waters  are  able  to  maneuver  to  a  suit¬ 
able  settlement  site  via  their  well-documented  selec¬ 
tive  tidal-stream  transport  mechanism  (reviewed  in 
Forward  et  al.  2003). 

We  divided  the  model  domain  into  14  estuaries  as 
shown  in  Fig.  1  based  upon  the  locations  where  juve¬ 
nile  and  adult  crabs  are  known  to  live.  These  estuar¬ 
ies  were  divided  based  upon  the  major  topographic 
feahues  in  the  region.  The  spawning  habitat  was 
divided  among  estuaries  according  to  geographic 
distance  as  described  above.  Each  larva  that  met  the 
other  criteria  for  settlement  was  deemed  to  settle  in 
the  nearest  estuary  as  defined  by  the  boundaries  in 
Fig.  1.  Neither  larval  development  nor  mortality  are 
incorporated  in  the  IBM;  however,  any  larvae  that  fail 
to  settle  in  an  estuary  in  60  d  are  assumed  to  either  be 
carried  by  alongshore  currents  over  the  western  or 
eastern  boundaries  of  the  model,  entrained  in  off¬ 


shore  currents  and  carried  to  the  deep  sea  off  the 
southern  boundary,  or  trapped  by  local  currents  off¬ 
shore  but  within  the  model  boundaries. 

Basic  dispersal  metrics  measured  or  tracked  inclu¬ 
ded  total  proportion  of  larvae  that  successfully  set¬ 
tled  into  estuaries,  fates  of  those  larvae  that  failed  to 
settle,  and  the  dispersal  distance  for  each  larva. 
Absolute  expected  larvae  dispersal  distance  from 
each  spawning  location  was  calculated  by  averaging 
the  migration  distances  of  successfully  settled  larvae 
in  aU  model  runs.  We  visually  examined  the  larval 
transport  trajectories  and  described  how  these  pat¬ 
terns  relate  to  published  descriptions  of  the  circula¬ 
tion  patterns  in  the  nGOM. 

Metapopulation  structure 

We  populated  24  potential  connectivity  matrices  P, 
whose  elements  ptj  are  the  proportion  (probability)  of 
larvae  from  source  estuary  i  that  settled  in  estuary  j 
for  each  simulated  spawning  date  f.  The  elements  of 
the  overall  potential  connectivity  matrix  P  are  the 
proportions  of  all  larvae  averaged  across  all  24 
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spawning  dates.  From  the  connectivity  matrices,  we 
also  populated  an  oceanographic  distance  matrix  D. 
The  elements  djj  are  given  by: 


where  s  is  the  step  discount  factor  that  accounts 
for  post-settlement  survival  and  recruitment.  The 
oceanographic  distance  dij  is  short  when  the  proba¬ 
bility  of  moving  from  estuary  i  to  estuary  j  is  high. 
The  probability  that  an  individual  vnU  follow  an  n 
step  pathway  through  multiple  estuaries  may  be 
computed  as  the  product  of  the  probabihtles  of  fol¬ 
lowing  each  step,  discounted  by  s”.  Alternatively,  the 
distance  of  this  pathway  may  be  computed  directly 
by  summing  the  distance  of  each  step.  For  example, 
a  pathway  from  estuary  i  to  j  to  k  would  have  proba¬ 
bility  ^PijPjk  and  distance  d,y  +  dy*.  Using  this  formu¬ 
lation,  it  is  possible  that  a  multiple  step  pathway  is 
shorter  and  more  probable  than  a  direct  pathway 
between  any  pair  of  estuaries,  Betweeimess  central¬ 
ity  for  a  given  estuary  is  the  proportion  of  minimum 
length  paths  in  the  network  that  pass  through  that 
estuary  (Freeman  1977,  1978).  Estuaries  with  a  large 
betweenness  centrality  are  important  for  maintaining 
connectivity  throughout  the  network.  We  measured 
betweenness  centrahty  for  each  estuary  from  the 
graph  generated  from  the  distance  matrix,  D,  using 
the  R  statistical  package  (R  Core  Team  2014)  with  the 
igraph  library  (Csardi  &  Nepusz  2006).  We  arbitrarily 
set  the  value  of  s  to  0.75. 

To  quantify  the  overall  level  of  coimectivity  in  the 
network  over  time,  we  computed  the  average  source 
strength,  average  vertex  degree,  and  Shannon  index 
for  each  vertex  from  each  potential  connectivity 
matrix  Pj.  The  source  strength  is  the  proportion  of 
particles  released  from  an  estuary  that  successfully 
settled.  This  metric  is  corrunonly  reported  in  ecologi¬ 
cal  network  studies.  Although  it  is  an  informative 
metric  for  assessing  larval  success  rates,  it  ignores 
the  number  of  possible  destinations  and  evenness 
in  connectivity  across  these  cormections.  Because 
mathematical  models  have  shown  that  cormectivity 
network  structure  is  an  important  driver  of  popula¬ 
tion  dynamics  (Kininmonth  et  al.  2010b),  we  use  2 
additional  metrics  of  population  cormectivity  along¬ 
side  source  strength.  The  vertex  degree  for  an  estu¬ 
ary  at  a  given  release  time  is  the  number  of  destina¬ 
tion  estuaries  that  received  settlers  from  that  estuary. 
Averaging  the  vertex  degree  across  all  estuaries  in 
the  network  provides  an  estimate  of  how  many  estu¬ 
aries  are  likely  to  be  impacted  by  a  perturbation  to  a 


single  estuary  within  1  generation.  The  Sharmon  index 
[Hi)  quantifies  both  the  number  and  strength  of  the 
outgoing  cormections  from  an  estuary  and  is  given  by; 

Hi  =  -J,Pii'in{pij)  (2) 

J 

The  Shannon  index  provides  an  estimate  of  the  sta¬ 
bility  of  the  metapopulation,  and  high  values  indicate 
a  more  diverse  suite  of  cormections  and  more  stable 
metapopulation.  We  averaged  each  of  these  3  metrics 
across  aU  14  estuaries  for  each  release  time  to  com¬ 
pute  a  measure  of  the  overall  connectivity  of  the 
population.  To  analyze  temporal  trends  in  connec¬ 
tivity  over  the  season,  we  performed  a  simple  linear 
regression  of  each  statistic  with  spawning  date. 

Deepwater  Horizon  oil  spill 

We  used  data  from  NOAA  Satellite  Derived  Sur¬ 
face  Oil  Analysis  Products  for  Deepwater  Horizon 
(www.ssd.noaa.gov/PS/MPS/deepwater.html)  to  de¬ 
termine  the  dciily  location  and  extent  of  oil  slicks, 
When  data  were  unavailable  for  a  specific  day,  we 
used  data  from  the  previous  day.  When  2  consecutive 
days  had  no  data  available,  we  used  data  from  the 
previous  day  for  the  first  day  and  data  from  the  fol¬ 
lowing  day  for  the  second  day.  If  a  larva  encountered 
an  oil  slick  at  any  point  during  its  trajectory,  we  des¬ 
ignated  that  larva  as  'oiled'.  We  quantified  the  poten¬ 
tial  decrease  in  connectivity  if  'oUed'  particles  were 
removed  from  the  simulation.  We  computed  the  total 
number  of  successfully  settling  larvae  that  were 
potentially  exposed  to  oU,  as  weU  as  the  location  of 
settlement  for  these  larvae.  From  these,  we  calcu¬ 
lated  the  total  proportion  of  settlers  to  each  estuary 
and  from  each  release  time  that  were  potentially 
exposed  to  oil.  Under  the  worst-case  scenario  that  ex¬ 
posure  to  oil  is  always  fatal  to  larvae,  we  reassessed 
the  spatial  structure  of  population  connectivity  in  the 
region.  Although  this  assumption  is  unlikely  to  be 
realistic,  it  provides  a  worst-case  bound  on  the  im¬ 
pacts  of  the  Deepwater  Horizon  oil  spill  on  blue  crabs 
during  the  dispersal  phase. 

Model  sensitivity  analyses 

To  confirm  that  our  results  are  robust  to  the 
assmnptions  made  in  the  model,  we  ran  a  number  of 
sensitivity  analyses.  Our  primary  sensitivity  analysis 
was  to  the  spawning  schedule.  There  is  some  evi¬ 
dence  for  limar  periodicity  in  blue  crab  spawning 
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(Ziegler  2002),  so  we  examined  the  scenario  in  which 
blue  crab  spawning  events  coincide  with  the  full  and 
new  moons.  In  our  analysis,  we  assume  that  the  max¬ 
imum  distance  migrated  by  ovigerous  female  crabs  is 
100  km.  There  is  no  strong  consensus  in  the  literature 
as  to  the  distance  migrated  by  female  crabs,  so  we 
tested  the  effects  of  doubling  and  halving  this  dis¬ 
tance  on  our  results. 


RESULTS 

Over  aU  time  periods,  7.7  %  of  the  simulated  larvae 
settled  within  the  region.  The  remaining  92.3  %  were 
removed  from  the  model  because  they  were  pre¬ 
dicted  to  move  out  of  the  model  domain  or  remain 
offshore  until  the  megalopae  were  no  longer  compe¬ 
tent  to  settle.  Of  the  total  48  million  larvae  simulated, 
34.8  million  (72.5%)  were  lost  to  the  edges  of  the 
model  domain  boundciries,  with  7.7  million  (22.2%) 
of  these  lost  over  the  southern  boundary,  1.7  million 
(4.9%)  over  the  eastern  boundary,  and  25.4  million 
(72.9%)  over  the  western  boundary.  This  pattern  of 
westbound  transport  being  dominant  was  also  exhib¬ 
ited  by  the  particles  that  successfully  settled;  75.7% 
of  the  successful  settlers  returned  to  their  origin 
estuary  or  one  to  the  west. 

The  trajectories  followed  by  the  larvae  were  in 
accordance  with  previous  descriptions  of  the  sea¬ 
sonal  circulation  of  the  northern  Gulf  of  Mexico.  A 
cyclonic  gyre  dominates  the  drcuiation  of  the 
Louisiana-Texas  (LATEX)  shelf  during  the  fall,  win¬ 
ter,  and  spring  (Smith  &  Jacobs  2005).  Most  larvae 
released  west  of  the  Mississippi  River  delta  during 
the  spring  and  fall  releases  were  entrained  by  the 
strong,  westward,  coastal  component  of  this  gyre  and 
transported  past  the  western  boundary  of  NGOM- 
NFS.  One  exception  to  this  pattern  is  particles 
released  immediately  southwest  of  the  delta  itself, 
many  of  which  were  entrained  by  a  small,  cyclonic 
eddy.  During  the  summer,  flow  on  the  LATEX  shelf 
is  predominately  eastward,  and  many  larvae  were 
transported  towards  the  delta,  and  then  offshore.  On 
the  eastern  side  of  the  delta,  most  individuals  re¬ 
leased  during  the  spring  were  transported  west 
towards  the  delta,  and  then  offshore.  This  pattern 
changed  in  the  summer  to  be  dominated  by  eastward 
alongshore  transport,  then  western  alongshore  trans¬ 
port  during  the  fall.  The  logio-transformed  expected 
distance  for  larvae  to  be  transported  from  the  spawn¬ 
ing  site  to  the  settlement  location  decreased  over  the 
course  of  the  year  (slope  =  -3.9  km  d"‘,  adjusted  = 
0.642,  p<  0.001). 


The  model  predicts  a  high  degree  of  local  retention 
of  blue  crabs  in  the  nGOM.  Overall,  28.5%  of  larvae 
that  successfully  settled  returned  to  their  parent  estu¬ 
ary  and  37.5%  settled  in  an  adjacent  estuary.  The 
high  local  retention  can  be  seen  by  the  strong  diago¬ 
nal  pattern  in  Fig.  2.  It  is  also  clear  from  Fig.  2  that 
the  MRD  acts  as  a  barrier  to  dispersal,  as  only  4  %  of 
the  settling  larvae  crossed  the  mouth  of  the  Missis¬ 
sippi  River.  Many  of  the  larvae  that  crossed  the  Mis¬ 
sissippi  River  were  swept  southward  offshore  and  did 
not  successfully  settle. 

Metapopulation  structure 

Total  potential  comiectivity  varies  across  estuaries 
and  over  time.  The  estuaries  adjacent  to  the  Missis¬ 
sippi  River  delta  have  the  highest  values  of  between¬ 
ness  connectivity,  indicating  that  they  are  most  im¬ 
portant  to  maintaining  the  connectedness  of  the 
metapopulation  (Fig.  3).  The  overall  connectivity,  as 
measured  by  the  Shaimon  index  and  vertex  out- 
degree  declines  over  the  course  of  the  spawning  sea¬ 
son.  Average  vertex  degree  drops  by  approximately 
0.029  links  per  day  (p  <  0.001)  and  the  Shannon  index 
by  0.004  per  day  (p  <  0.001;  Fig.  4).  A  non-significant 
(p  5  0.13)  increase  was  seen  in  source  strength. 

Deepwater  Horizon  oil  spill 

Of  the  3.7  million  successfully  settling  larvae, 
1.4  million  (38.1%)  were  potentially  exposed  to  oil 
during  the  dispersal  phase.  The  proportion  of  suc- 
cessfiilly  settling  larvae  that  were  potentially  ex¬ 
posed  to  oil  varied  widely,  from  0  %  when  there  was 
no  overlap  to  93%  for  larvae  released  on  18  June 
2010,  when  the  overlap  was  greatest.  The  majority  of 
the  potentially  oiled,  successfully  settling  larvae 
were  released  during  June,  with  13.8%  on  12  Jime 
2010  and  20,9%  on  26  June  2010.  The  overlap  of  the 
oil  slick  with  larval  locations  was  much  greater  on  the 
eastern  side  of  the  MRD  than  the  western  side,  so  it  is 
not  surprising  that  79.5%  of  the  settled  larvae  that 
were  potentially  exposed  to  surface  oil  were  pre¬ 
dicted  to  settle  in  estuaries  east  of  the  MRD  and  only 
20.5%  in  estuaries  west  of  the  MRD  (Fig.  5).  On 
26  Jime  2010,  96.3  %  of  the  successM  settlers  east  of 
the  MRD  were  potentially  exposed  to  oil. 

Larvae  were  most  likely  to  be  exposed  to  oil  at  a 
young  age  and  for  a  short  duration:  63.2%  of  the 
potentially  exposed  larvae  would  have  been  first 
exposed  at  less  than  15  d  of  age,  and  28.6%  were 


Source  Source 

(Number  settled  in  1000s)  (Number  settled  in  1000s) 
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potentially  exposed  for  only  a  single 
day.  The  distribution  of  exposure 
durations  was  approximately  expo¬ 
nentially  distributed,  with  a  decay 
rate  of  0,34  d“^. 

Assuming  that  exposure  to  oil  is 
fatal  to  blue  crab  larvae  reduces  the 
proportion  of  larvae  that  success¬ 
fully  settle  by  38.1%,  but  does  not 
substantially  influence  the  patterns 
of  cormectivity  (Fig.  2).  Overall,  the 
proportion  of  successfully  settling 
larvae  from  each  estuary  was  re¬ 
duced  by  31.8  ±  15.3%.  Larvae  re¬ 
leased  from  the  2  estuaries  border¬ 
ing  the  MRD  were  particularly 
heavily  impacted  by  the  oil  spiU, 
■with  47.1%  of  the  larvae  leaving 
Adams  Bay  and  51.6%  of  the  larvae 
leaving  Spanish  Pass  being  poten¬ 
tially  exposed  to  oil.  As  recipients, 
Barataria  Bay,  Perdido  Bay,  and 
Chocawhatchee  Bay  were  most 
heavily  impacted  with  78.7%, 
69.7%,  and  67.1%  of  their  settlers 
potentially  exposed  to  oil.  These 
sites  were  particularly  impacted  be¬ 
cause  they  already  received  few  set¬ 
tlers  and  much  of  their  settlement 
occurred  during  the  DWH  oil  spill. 
Other  sites  that  consistently  had 
high  settlement  throughout  the  year 
(e.g.  Chandeleur  Sound)  were  less 
impacted  by  the  spUl.  If  oil  exposure 
were  fatal,  Chandeleur  Sound 
would  be  more  important  for  main¬ 
taining  connectivity  (Fig.  3). 

Because  the  oil  slicks  from  DWH 
were  present  during  late  spring  and 
early  summer,  fatality  from  oil  woiild 
have  dampened  the  trend  that  con¬ 
nectivity  declines  during  the  year 
(Fig.  4).  Excluding  potentially  oiled 
larvae,  the  Shannon  index  and  ver¬ 
tex  degree  declines  over  the  course 
of  the  spawning  season,  but  not  by 
as  much  as  with  those  larvae  in¬ 
cluded,  The  average  vertex  degree 
drops  by  approximately  0.016  links 
per  day  (p  <  0.001)  and  the  Shannon 
index  by  0.003  per  day  (p  <  0,001). 
Although  previously  non-significant, 
the  increase  in  source  strength  dur- 
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section  of  the  text.  This  metric  shows  the  importance  of  each  estuary  to  maintaining  connectivity.  (A)  All  larvae;  (B)  excludes 
larvae  that  were  potentially  exposed  to  oU.  In  both  cases,  Chandeleur  Sound  is  the  most  important  estuary  for  maintaining 
connectivity.  When  particles  that  are  potentially  exposed  to  oil  are  removed,  the  importance  of  Chandeleur  Sotmd  increases 


Fig.  4.  Average  vertex  degree  and  average  Shannon  index  for  each  simulation  versus  spawning  date  measuring  the  overedl 
connectivity  of  the  blue  crab  population  in  the  network  of  estuaries.  Error  bars  are  ±1  SE.  (A)  All  larvae;  (B)  excludes  larvae 
that  were  potentially  exposed  to  oil.  The  significant  trend  lines  for  average  vertex  degree  and  Shannon  index  are  plotted 
in  both  panels.  The  trend  for  source  strength  is  included  only  in  (B)  because  it  was  not  significant  when  oiled  particles 

were  included 


ing  the  year  is  significant  when  potentially  oiled  lar¬ 
vae  are  excluded  (p  <  0.005).  However,  the  rate  of 
increase  is  small,  5.4  x  10"*  d"V  T^ese  results  support 
our  conclusions  from  examining  the  coimectivity 
matrices  that  the  magnitude,  but  not  patterns,  of  con¬ 
nectivity  were  influenced  by  the  DWH  oil  spiU. 


Sensitivity  analyses 

Our  sensitivity  analysis  to  the  frequency  of  larval 
release  revealed  that  our  model  is  not  sensitive  to  this 
configuration.  Assmning  that  blue  crabs  spawn  only 
during  the  full  and  new  moons  reduced  the  number 


Average  Shannon  index  and  source  strength 
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Fig.  5.  Successfully  settling  blue  crab  larvae  from  each  release  time  divided  according  to  the  proportion  that  arrives  at  each 
estuary.  Grey  bars  indicate  larvae  that  were  not  exposed  to  oil,  and  black  bars  indicate  larvae  that  were  potentially  exposed 
to  oil  at  any  time  during  the  dispersed  phase.  The  dashed  line  is  the  location  of  the  Mississippi  River  delta.  Estuaries  are  listed 
from  west  to  east  (left  to  right).  WCB:  West  Cote  Blanche  Bay;  ATB:  Atchafalaya  Bayj  CB:  Caillou  Bay:  TB:  Terrebonne  Bay; 
BB:  Barataria  Bay;  AB:  Adams  Bay:  SP;  Spanish  Pass:  CS:  Chandeleur  Sound;  G:  Gulfport:  P;  Pascagoula;  MB:  Mobile  Bay;  PB: 

Perdido  Bay;  PSB:  Pensacola  Bay;  CHB:  Chocawhatdiee  Bay 
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of  successful  larvae,  but  not  the  probability  of  suc¬ 
cess.  In  our  original  simulation,  7,7  %  of  the  released 
larvae  successfully  settled.  Examining  the  new  and 
full  moon  events  only,  this  proportion  was  8.2%.  The 
overall  patterns  that  the  MRD  is  a  barrier  to  dispersal 
(Fig.  SI  in  the  Supplement  at  www-int-res.com/ 
articles/suppl/m527pl43_supp.pdf)  and  that  the  es¬ 
tuaries  on  either  side  of  the  MRD  are  most  important 
for  maintaining  connectivity  (Fig.  S2  in  the  Supple¬ 
ment)  remained.  As  before,  we  saw  a  significant 
decline  in  connectivity  as  measured  by  the  Shannon 
Index  (p  <  0.001)  and  average  vertex  degree  (p  S 
0.001). 

Reducing  the  distance  migrated  by  ovigerous 
females  to  50  km  and  doubling  it  to  200  km  does  not 
change  our  results  (Fig,  S3  in  the  Supplement),  The 
MRD  remains  a  barrier  to  dispersal,  and  the  trends 
in  Shannon  index,  average  vertex  degree,  and 
source  strength  are  unchanged  in  all  configurations. 
In  our  original  parameterization,  average  vertex  de¬ 
gree  drops  by  -0.029  links  per  day  (p  S  0.001)  and 
the  Shannon  index  by  0.004  per  day  (p  <  0.001),  Our 
results  here  are  highly  similar.  Under  the  50  km 
parameterization,  average  vertex  degree  drops  by 
-0.028  links  per  day  (p  <  0.001)  and  the  Shannon 
index  by  0.003  per  day  (p  S  0.001).  Under  the 
200  km  parameterization,  average  vertex  degree 
drops  by  -0.029  links  per  day  (p  S  0.001)  and  the 
Shannon  index  by  0.004  per  day  (p  <  0.001).  A  non¬ 
significant  increase  was  seen  in  source  strength  in 
our  original  configuration  (p  <  0.13),  the  50  km  con¬ 
figuration  (p  <  0.12),  and  the  200  km  configuration 
(p  S  0;13).  Overall,  our  results  are  robust  to  the 
assumptions  that  we  have  made. 

DISCUSSION 

The  Mississippi  River  and  its  delta  have  a  strong 
effect  on  the  circulation  in  the  nGOM,  Strong  cur¬ 
rents  often  occur  near  the  Mississippi  River  along 
the  river  plumes  that  advect  larvae  away  from  the 
river  mouth,  and  alongshore  currents  form  on  either 
side  of  the  delta,  driven  by  the  wind  and,  to  a  lesser 
extent,  by  the  tides.  Hence,  the  MRD  forms  a  natu¬ 
ral  barrier  to  dispersal  (Fig.  2).  This  has  a  large 
effect  on  the  potential  population  connectivity  of 
blue  crabs  within  the  nGOM.  Our  results  suggest 
that  the  population  east  of  the  delta  might  be  eco¬ 
logically  distinct  from  the  population  west  of  the 
delta.  No  evidence  exists  for  genetic  differentiation 
at  the  MRD  for  blue  crabs,  although  a  review  of 
genetic  studies  in  the  Gulf  does  show  that  other 


species  differentiate  genetically  at  the  MRD  (Neigel 
2009).  However,  differentiation  in  benthic  species 
could  be  due  to  sediment  differences  east  and  west 
of  the  MRD  or  other  factors,  rather  than  larval  trans¬ 
port  (Neigel  2009).  Furthermore,  the  small  cross¬ 
delta  movements  predicted  by  our  model  may  be 
large  enough  to  homogenize  genetic  differences 
(Slatkin  1987),  and  we  do  not  necessarily  expect  to 
find  genetically  differentiated  populations.  On  eco¬ 
logical  time  scales,  the  larval  connectivity  patterns 
found  here  support  managing  each  side  of  the  MRD 
separately. 

We  speculate  that  the  decline  in  population  con¬ 
nectivity  during  the  summer  is  primarily  driven  by 
the  seasonal  circulation  patterns  of  the  LATEX 
shelf.  Whereas  both  average  vertex  degree  and  the 
Shannon  index  indicated  a  significant  decline  in 
population  connectivity  during  the  year,  the  source 
strength  did  not  show  a  significant  trend.  Theoreti¬ 
cal  studies  have  found  that  network  structure  is  an 
important  driver  of  metapopulation  dynamics  (Hast¬ 
ings  &  Botsford  2006,  Kininmonth  et  al.  2010b),  so 
we  recommend  that  future  studies  of  population 
connectivity  apply  a  suite  of  connectivity  metrics. 
The  trajectories  followed  by  the  simulated  larvae 
are  similar  to  published  circulation  patterns  over  the 
LATEX  shelf.  Seasonal  variability  in  potential  con¬ 
nectivity  determines  how  timing  of  spawning  peaks 
■within  the  year  affects  dispersal  distances  and  real¬ 
ized  population  cormectivity.  Because  there  is  some 
interannual  variability  in  the  circulation  patterns, 
particularly  on  the  eastern  side  of  the  MRD  (Smith  & 
Jacobs  2005),  and  because  stochasticity  in  realized 
connecti'vity  lowers  metapopulation  growth  rates 
(Watson  et  al.  2012),  it  would  be  interesting  and 
useful  to  quantify  the  interannual  variability  in  con¬ 
nectivity  in  this  region. 

Our  model  can  only  be  used  to  predict  potential 
connectivity  as  it  uses  assumptions  of  constant  larval 
supply  and  does  not  accoimt  for  larval  mortality  dur¬ 
ing  transport  or  mortality  at  later  stages.  To  predict 
realized  connectivity  would  require  knowledge  of 
how  the  larval  supply  varies  over  space  and  time  as 
well  as  more  accurate  information  on  spawning  loca¬ 
tions,  larval  mortality,  finer-scale  nearshore  hydrody¬ 
namics,  and  post-larval  beha'viors  that  affect  settle¬ 
ment.  Such  information  does  not  exist  for  blue  crabs 
in  the  nGOM.  Since  our  dispersal  model  does  not 
include  such  factors,  we  do  not  expect  resulting  set¬ 
tlement  patterns  predicted  by  our  model  to  match 
data  from  megalopal  (post-larval)  collections  (e.g. 
Rabalais  et  al.  1995,  Perry  et  al.  2003,  E.  K.  Grey  et  al. 
unpubl.  data). 
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Use  of  graph-theoretic  metrics  improves  the  analy¬ 
sis  of  results  from  Lagrangian  particle-tracking  mod¬ 
els.  It  can  be  argued  that  not  much  additional  infor¬ 
mation  is  gained  by  calculating  the  graph-theoretic 
metrics  for  such  a  small  network  (only  14  vertices), 
especially  when  the  vertices  have  an  essentially 
linear  geographic  distribution  along  the  coast.  How¬ 
ever,  we  believe  that  the  use  of  graph  metrics  in  a 
small  network  demonstrates  the  effectiveness  and 
utility  of  the  approach  and  of  the  particular  metrics. 
For  example,  we  expect  that  the  estuaries  immedi¬ 
ately  adjacent  to  the  dispersal  barrier  imposed  by  the 
MRD  win  be  important  stepping  stones  since  the  low 
percentage  of  larvae  that  travel  across  the  MRD 
largely  originate  from  these  estuaries.  Betweenness 
centrality  is  indeed  highest  for  those  estuaries 
(Fig.  3),  demonstrating  that  betweenness  centrahty  is 
an  effective  tool  for  measuring  the  importance  of 
individual  vertices,  and  therefore  can  be  used  in 
larger  networks  where  the  geographic  arrangement 
of  vertices  is  more  complex  (e.g.  Treml  et  al.  2008, 
Kininmonth  et  al.  2010a,  Treml  &  Halpin  2012).  Simi¬ 
larly,  average  vertex  degree  and  the  Shannon  index 
are  both  effective  metrics  for  measuring  the  overall 
connectivity  of  networks  (Shannon  1948,  MacArthur 
1955). 

The  DWH  oil  spill  overlapped  with  predicted  blue 
crab  spavraing  locations  in  late  May  and  June  2010. 
Our  model  predicts  that  larvae  potentially  exposed  to 
oil  are  most  likely  to  settle  east  of  the  MRD,  with 
peak  settlement  in  late  June.  No  catastrophic  reduc¬ 
tion  in  settlement  is  observed  in  this  or  any  other 
region  during  this  time  period  (E.  K.  Grey  et  al. 
unpubl.  data).  However,  little  is  known  about  the  sub- 
lethal  effects  of  oil  exposure  in  crabs.  If  effects  mani¬ 
fest  later  in  life,  this  model  predicts  that  we  would 
observe  such  effects  more  commonly  in  eastern  nGOM 
estuaries. 

Overall,  we  find  that  the  MRD  serves  as  a  barrier  to 
larval  dispersal  of  blue  crabs  in  the  nGOM,  and  that 
the  overall  level  of  population  connectivity  declines 
during  the  spawning  season.  Our  results  are  robust  to 
the  assmnptions  that  we  have  made  with  regards  to 
the  distance  that  ovigerous  female  crabs  migrate  and 
the  frequency  of  spawning.  Blue  crabs  that  were  po¬ 
tentially  exposed  to  oil  from  the  DWH  oil  spill  are  pre¬ 
dicted  to  have  settled  on  the  eastern  side  of  the  MRD. 
This  study  provides  a  valuable  first  estimate  of  the  po¬ 
tential  impacts  of  the  DWH  oil  spill  on  larval  transport. 
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